Nuclear/Circumnuclear Starbursts and Active Galactic Nuclei 
Mass Accretion in Seyfert Galaxies 



Yasuyuki Watabe 1 ' 2 '*, Nozomu Kawakatu 3 and Masatoshi Imanishi 3 



1 Center for Computational Sciences, University of Tsukuba, Ten-nodai, 1-1-1 Tsukuba, 

Ibaraki 305-8577, Japan; watabe@ccs.tsukuba.ac.jp 
2 INAF-Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy 
3 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181- 

Japan 



ABSTRACT 

We investigated the correlation between nuclear /circumnuclear starbursts 
around the active galactic nuclei (AGNs) and the AGN activities for 43 Seyfert 
galaxies in the CfA and 12 /zm samples. We found that circumnuclear starburst 
luminosity as well as nuclear starburst luminosity are positively correlated with 
AGN luminosity Moreover, nuclear starburst luminosity is more strongly cor- 
related with the AGN luminosity normalized with AGN Eddington luminosity 
than is circumnuclear starburst luminosity This implies that starbursts nearer 
the AGN could have a greater effect on AGN mass accretion. We also discuss 
these results from the viewpoint of the radiation effects from starbursts and se- 
quential starbursts. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert - 
galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 



Since the discovery of active galactic nuclei (AGNs) the physical mechanism of AGN fuel 
ing remains unresolved. Various fueling mechanisms ha ve been considered thus far, for exam- 



ple, t i dal torque driven by majo r /minor galaxy merger (IHernouist 



19911 : iMihos fc Hernquistl 11990 : iTaniguchil Il999t ISaitoh fc Wada 



nonaxisymmetric gravitational potential due to a stellar bar (I No guchil 1 1 9 8 8l : IShlosman et al. 



989l: iBarnes fc Hernquist 



2004) , tidal torque from 
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1990l ; Barnes fa Hernquistlll992l ; iKnapen et alll995 1 Benedict et al 



shock or turbulence in the interstellar medium 



Wada fa Normanlll999l . l200ll ; iMacieiewski et al. 



Fukuda et al 



2002 



and radiation drag from a starburst (jUmemura et al 
19981 : lohsuga et al.lll99gh . 



Wada et al 



1997 



1996 



Fukuda et al. 



1998: IMontenegro et al 



2002: Qnodera et al. 



Fukue et al.l 1997: Umemura et al. 



200 



3), 



1999; 



2004), 



Numerous observations have gradually clarified circumnuclear starburst events (several 



ten exhibit patchy and ringlike structures 


(Pogge! 


198 


9; 


Wilson et al. 199 


1: Forbes et al. 


1994; 


Marconi et al.l 


1994: Mauder et al. 1994: ] 


3uta et al. 


19951: Barth et al. ] 


L995|; 


Leitherer et al. 


1996; 


Maoz et al. 


1996: Storchi-Bergmann et al.l 


1996; 


Elmouttic ct al. 


1998 


Knapen et al. 


2002; 


Knapen 


20051). In addition, hidden nuclear starbursts (within a few hundred par- 



2003 



sees from the center) are found both in Seyfert 1 and 2 galaxies (llmanishil 12002 
Rodriguez- Ardila fa Viegasl 120031) and the nuclear s tarburst luminosity is positively corre 
lated with the AGN power (llmanishi fa Wadal 120041 ) . 



These starburst events could strongly influence the s tructure and the dynamics of gas 



through the energy input of mu l tiple supernova explosion s ( Shapiro fa Fieldfl976 : Tomisaka fa Ikeuchi 

nl 2002 ) and the effects of strong radiation 



1986 



Norman fa Ikeuchi 



Umemura et al 



1997 



Wata 



je fa Umemura 



1989 



1998 



2005 



Wada fa Normar 



Qhsuga fa Umemura 



19991 . 120011 : iKawakatu fa Umemural 12002 



Thompson et al.l 120051 ) . Especially, for the AGN mass accre- 



Wada & Norman J2002l )~ shows the possibility of AGN fueling from the effect of the 



tion, 

chaotic disturbance by the nuclear starburst us ing three-dimensional hydrodynamic sim- 
ulation. Moreover, lUmemura et al.l (119971 . Il998l ) have shown that the radiation effects of 
circumnuclear starburst could cause AGN fueling. 

These theoretical and observational results suggest that some connections exist between 
the nuclear /circumnuclear starburst and AGN activi ty. Although a relationship between 
the n uclear starburst and the AGN has been reported (llmanishi fa Wadall2004l : iDavies et al. 
20071 ). this remains uncertain regarding the circumnuclear starburst. Moreover, it is not clear 
whether a specific relationship exists between AGN mass accretion and the starburst at vari- 
ous scales. In this study, we investigated the relationship between the nuclear /circumnuclear 
starburst and AGN activity by considering both nuclear and the circumnuclear starbursts. 

To investigate nuclear /circumnuclear starbursts, we used the polycyclic aromatic hy- 
drocarbon (PAH) emission feature at 3.3 /im, 6.2 /im, 7.7 /im, and 11.3 /xm. The PAHs 
are excited by far-UV photons and strong PAH emission are seen in star-forming regions, 
whereas a pure AGN shows only a featureless spectrum with virtually no PAH emission. Also, 
since the dust extinction is much lower for these 3.3 — 11.3/zm PAH emission (< 0.06Ay; 
Rieke fa Lebofskylll985l : iLutz et al.lll996l ). we can quantitatively estimate modestly obscured 
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(Ay < 15 mag) star-formation activity by using these emission lines. Throughout the paper, 
we adopted H = 80 kms~ 1 Mpc~ 1 , Qm = 0.3, and = 0.7. 



SAMPLE DATA 



We studied 2 1 Seyfert 1 galaxi es and 22 Seyfert 2 galaxies in the CfA ( jHuchra fc Burg 



19921 ) and 12 /xm (IRush et al.lll993l ) samples selected based on their host-galaxy magnitudes 



and IRAS 12 fim fluxes, respectively. These samples are not expected to b e biased toward or 



agains t the presence of nuclear starbursts. Our sample was selected from llmanishi fc Wada 



(120041 ) data set with both the 3.3 /im PAH luminosity and the nuclear N-band luminosity. 
In their sample, Seyfert galaxies at z — 0.008 — 0.035 were selected to investigate the nuclear 
starburst by ground-based spectroscopy using a l"-2" slit, where 1" corresponds to a physical 
scale of 150 (z = 0.008) to 650 pc (z = 0.035). Also, to be best observable from Mauna Kea, 
Hawaii, the declinations of Seyfert galaxies are limited to being larger than —30°. Owing to 
the telescope limit of the IRTF 3 m telescope, a restriction of declination of less than 68° is 
also applied. For their sample, there are no obvious biases. 

Among 43 Seyfert galaxies, we estimated the supermassive black hole (SMBH) mass 
to investigate the mass accretion rate normalized with the AGN Eddington mass accretion 
rate for 25 objects. To these objects, we estimated the 6.2 /im, 7.7 /im, and 11.3 /im PAH 
emission for 13 objects further. Also, we collected the X-ray luminosity from the data in 
the literature for 15 objects. Since our sample included some upper limit data, we used 
statistical techniques applicable to a sample where upper limit data are present. 

We focused on both nuclear starbursts, which exist within a few hundred parsecs from 
the center, and circumnuclear starbursts in the entire host galaxy region. Figure [1] shows the 
schematic view of the nuclear /circumnuclear starbursts. We can consider the 3.3 /im PAH 
emission, which were derived through ground-based spectroscopy using the narrow slit, as 
a nuclear starburst. Also, since a circumnuclear starburst is typically much higher than a 
nuclear starburst (see section 4.2), the 6.2 /im, 7.7 /im, and 11.3 /tm PAH emission, which 
were observed in the whole host galaxy region, represent a circumnuclear starburst. In detail, 
we described in section 3.1 and 3.2. 
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ESTIMATION OF PHYSICAL PROPERTIES 



3.1. Nuclear Starburst 



P AH molecules near the AGN can be destroyed by strong X-ray radiation from the 
AGN ( )Voitlll992l ; ISiebenmorgen et al.ll2004l ). However, if the PAH molecules are sufficiently 
shielded by a substantial column density of X-ray-absorbing gas, PAHs can survive even if 
they exist near the AGN. In fact, ample gas is believed to exist that could be related to 
the obscuring material (the dusty torus in the context of the AGN unified model) around 
the AGN. In addition, the gravitational stability parameter Q ( )Toomrelll964l ) decreases with 
the radius of th e obscuring material, whose mass is notably smaller than that of the SMBH 
( llmanishil 120031 ) . The gravitational collapse of molecular gas can therefore occur more easily 
in the outer part of the obscuring material. Thus, in the outer part of the obscuring material, 
a starburst could exist near the AGN and PAHs could survive. In fact, from the estimation 
of the surface brightness val ues of the 3.3 //m PAH emission (jlmanishi fc Alonso-Herrero 
20041 ; llmanishi &: Wadall2004l ). which were derived through ground-based spectroscopy using 
the l"-2" slit, this emission should come from nuclear starbursts in the gas-rich region, which 
may be the obscuring material around AGN. 

Therefore, we can use their 3.3 /xm PAH emission data, which were observed by ground- 
based spectroscopy using the narrow slit, as the indicators of nuclear starbursts. Although 
the circumnuclear starburst, which often exhibit ring-like structures with a radius of kilo- 
parsecs order from the center, along the slit direction could be also included, the fraction 
of the circumnuclear starburst emission inside this thin slit is negligible compared to the 
whole. Thus, the 3.3 /xm PAH emission could be regarded as a practical probe for nuclear 
starbursts in both Seyfert 1 and Seyfert 2 galaxies. 



In llmanishil (120031 ) and llmanishi fe Wadal (120041 ). only upper limits for the nuclear 3.3 
jum PAH fluxes are available in more than half of the observed Seyfert galaxies. Additionally, 
the observed Seyfert galaxies do not comprise a complete sample. To increase the number 
and fraction of Seyfert galaxies with detectable nuclear 3.3 //m PAH emission, we newly 
observed eight sources, which are shown in Table HJ Among them, six sources (NGC 931, 
F03450+0055, NGC 262, NGC 513, MCG-3-58-7, and Mrk 993) showed poss ible signs of the 



3.3 /xm PAH emission i n our previously obtained infrared 3-4 /xm spectra (llmanishil 12003 
Imanishi fc Wadal 120041 ) . Therefore, we decided to re-observe these sources with longer 



exposure to check if the signs are real or not. The two sources (MCG-2-8-39 and 0152+06) 
are newly observed this time to increase the number of observed sources. We showed these 
observations in Appendix. 
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3.2. Circumnuclear Starburst 



To estimate the absolute magnitudes of circumnuclear starburst activity in Seyfert 
galaxies, we analyz ed archival infrared data obtained with the Infrared Spectr ograph (IRS) 
( IHouck et al.ll2004l ) onboard the Spitzer Space Telescope (IWerner et al.l I2004J ) through the 
program PID 3269 (PI = Gallimore). We focused on the infrared spectra of Short-Low 2 
(SL2; 5.2-7.7 /im) and 1 (SL1; 7.4-14.5 /im), to estimate the fluxes of PAH-emission features 
at A rcst = 6.2 /im, 7.7 //m, and 11.3 /im in the rest frame. 

We used the latest pipeline-processed data products at the time of our analysis (Sll-14, 
pbcd files). The observations were performed with the slit-scan mode, and the entire host 
galaxy regions of the Seyfert galaxies were covered. We summed all the signals and extracted 
the spectra after subtracting the background emission in a standard manner. Thus, the re- 
sulting spectra should reflect the emission from the entire regions of the Seyfert galaxies. 
Wavelength and flux calibrations for the SL data were made on the basis of the files of 
the Spitzer pipeline-processed data, termed "bO_wavsamp.tbl" and "bO_fluxcon.tbl" , respec- 
tively. For SL1 spectra, data at A b s > 14.5 /im in the observed frame are invalid (Infrared 
Spectrograph Data Handbook Version 1.0) and so were removed. For some faint sources, 
appropriate spectral binning was applied to reduce the scatter of data points and to achieve 
more reliable measurements of the PAH fluxes. 

The circumn uclear starburs t luminosity is typically much higher than the nuclear star- 
burst luminosity (llmanishill2003l ; llmanishi fc Wadall2004l . and see section 4.2). Thus, we can 
consider that the 6.2 /im, 7.7 /im, and 11.3 /im PAH emission, which were observed in the 
whole host galaxy regions, are good indicators of a circumnuclear starburst even if these 
PAH luminosity contain a small contribution from nuclear starbursts. 



3.3. AGN Luminosity 



For AGN activity, we used t he nuclear N-band (10.5/im) luminosity measured with a 
1.5" aperture (IGoriian et al.ll2004Tl since the contribution of starbursts for the nuclear N-band 
luminosity is small (IGorjian et al.ll2004l ). However, nuclear N-band luminosity represents the 
reemission of hot dust illuminated by the AGN. So it may contain some uncertainties due to 
the dust distribution and the covering factor of the dusty torus around the AGN. To examine 
whether the nuclear N-band luminosity represents AGN activity in the sample data, we also 
used absorption-corrected hard (2-10 keV) X-ray luminosity, which is regarded as the AGN 
intrinsic luminosity for the Compton-thin AGN by using the data in the literature (see 
references in Table |3]) . 
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3.4. Black Hole Mass 



Next, we estimated the SMBH mass, Mbh, to investigate the mass accretion rate nor- 
malized with the AGN Eddington mass accretion rate. To estimate SMBH ma sses of Seyfert 
1 galaxies, we used the following method (in detail, see lKawakatu et al.ll2007l ). The method 
for estimating the mass of a SMBH is based on the assumption that the motion of ionized 
gas clouds moving around the SMBH is dominated by the gravitational force and the clouds 
within the broad-line regi on (BLR) are viria lized (e.g.. iPeterson fc Wandellll999l ). Adopting 
an empirical relationship (IKaspi et al.ll2000l ) between the size of the BLR and the rest-frame 
optical continuum luminosity, AL,\(5000A) and reverberation mapping, we can obtain the 
following formula: 



Mi 



BH 



4.9+°; 4 x 10 6 



AL A (5100A) 
10 44 ergss- 1 



0.70±0.033 



^FWHM 



.10 3 kms 1 . 

We can use this equation for estimating the mass of SMBHs in Seyfert 1 galaxies. 



(1) 



We cannot use the above method for Seyfert 2 galaxies because no broad emission 
components exist. So we used the data of SMBH mass (IBian fc Gull2007l), which is esti mated 
by SMBH mass and the stellar velocity dispersion, cr*, relation (jTremaine et al.ll2002l ). 



Mi 



BH 



10 



8.13 



/ 0* 

V200kms 



4.02 



M„. 



(2) 



4. RESULTS 

4.1. Nuclear 3.3 pm PAH Emission 

Figure [2] represents the zoom-in spectra around the 3.3 /im PAH emission feature of the 
observed Seyfert galaxies. Thanks to the excellent observing conditions and longer exposure, 
the flux excesses at the expected wavelength of the 3.3 /im PAH emission feature at (1 + z) 
x 3.3 /im are now recognizable in the three sources (F03450+0055, NGC 513, and MCG-3- 
58-7), marked with "3.3 /im PAH" in Fig. El For these sources, the fluxes, luminosities, and 
rest-frame equivalent widths (EW 3 3 pah) of the 3.3 /im PAH emission were estimated using 
Gaussian fits. For the remaining five sources with no clear 3.3 /im PAH detection (marked 
with "3.3 /im PAH (?)" in Fig. [2]), we estimated the upper limits of the PAH strengths, by 
adopting the lowest pl ausible continuum levels and assuming the profile of type-1 sources in 



Tokunaga et al.l (119911 ) as a 3.3 /im PAH template. Table [2] summarizes the results. 



At A rest ~ 3.3 /im in the rest frame, PiS (3.30 /im) emission is present, superposed on 
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the 3.3 /im PAH emission. The relative contribution from this Pf<5 emission line is expected 
to be high in Seyfert 1 galaxies, because broad emission line components are unobscured. 
However, the equivalent width of the Pf<5 line is expected to be ~0.3 nm for a typical Seyfert 
1 galaxy (jlmanishi et al.ll2006l ). and even lower in a Seyfert 2 galaxy. This equivalent width 
value is much lower than the observed values or upper limits in Table [2j We thus ascribe 
the flux excess at A rcst = 3.3 /im mostly to the 3.3 /im PAH emission feature. 

Among the six previously observed Seyfert galaxies, the de tected 3.3 /xm PAH fluxes or 
their upper limits are lower than previously derived upper limits (Ilmanishill2003l ; llmanishi fc Wada 
2004h . except NGC 931. For NGC 931, our new upper limit of the 3.3 /im PAH flux is ~50% 
higher than our previous estimate (jlmanishi fc Wadall2004l ). possibly because our new esti- 
mate is a conservative one, derived based on the lowest plausible continuum level. The 3-4 
/im continuum flux levels a re also similar within ^50% between t he spectra in this paper 
and in our previous papers (Ilmanishill2003l ; llmanishi fc Wadal 120041 ) , except NGC 5 13. The 
flux of NGC 513 in this paper is more than 1 mag fainter than in llmanishil (120031 ). possi- 
bly because of narrower slit employed and/or the flux variation of the 3-4 /im continuum 
emission, originating in AGN- heated hot dust. 



4.2. 6.2 /im, 7.7 /im, and 11.3 /mi PAH Emission 



Figure [3] and H] represent 5.2-14.5 /im spectra of Seyfert 1 and 2 galaxies, respectively. 
Most of the Seyfert galaxies in Figure [3] and H] show clearly detectable PAH emission features 
at A rest = 6.2/im, 7.7/im, and 11.3/zm. 

To estimate the fluxes of the PAH emission features at A rest = 6.2 /im, 7.7 /im, and 11.3 
/im, we adopted a linear continuum determined from d ata points at the shor ter and longer 
wavelength sides of individual PAH features, as used by llmanishi et al.l (120071 ). The adopted 
continuum levels are shown as solid straight lines in Fig. [3] and HJ We fitted the PAH 
emission features with Gaussian profiles, which reproduced the observed data reasonably 
well. We treated the weak (lower than 3 a) or undetected PAH emission as the upper limit 
or non-flux, respectively. 

By using 3.3 /im, 6.2 /im, and 11.3 /im PAH luminosities, we can roughly estimate 
infrared luminosities of the nuclear starburst and total starburst. In Table HJ we showed 
infrared luminositie s of the nuclear starburst estima ted from the 3.3 /im PAH luminosity; 
-^ir,3.3 — -^3.3 x 10 3 ( Mouri et al. 1990 ; ImanishJ 2002 ) and the total starburst estimated from 
the 6.2 /im and 11.3 /im PAH luminosity; Z /TR.fi.2 = -^6.2 x 3 x 10 2 ( IPeeters et al.l 120041 ) and 
-^ir,ii.3 = £11.3 x 7 x 10 2 (jSoifer et al.l 120021 ). respectively. Infrared luminosities of the total 
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starbursts are several times or more as large as those of nuclear starbursts. This means that 
the circumnuclear starbursts dominate the total starburst luminosity. Thus, we can consider 
the 6.2 /im, 7.7 fim, and 11.3 /im PAH emission, which were observed in the whole host 
galaxy regions, as good indicators of a circumnuclear starburst. 



4.3. Luminosity Correlation between Nuclear/Circumnuclear Starburst and 

AGN 

In Figure we checked a correlation between the nuclear N-band luminosity and 



tics ( 


Isobe et al. 


1986) 


Feieelson & Nelson 


1985 



rovided in the A stronomy Survival Analysis package (ASURV; 



Isobe et al.l Il986l ) to both types of Seyfert galaxies. The prob- 
ability that a correlation is not present was 0.9 % (2.6 a). Since these hard X-ray luminosity 
are the absorption-corrected luminosity and estimated for the Compton-thin AGN, this cor- 
relation means that the nuclear N-band luminosity is good indicator of AGN activity for 
these Seyfert galaxies samples. 

In Figure [61 we plotted the nuclear N-band luminosity and the hard X-ray luminosity 
versus the 3.3 /im PAH emission luminosity. The probability that a correlation is not present 
was 0.7 % (2.7 a) and 5 .3 % (1.9 a) for the nucl ear N-band luminosity and the hard X-ray 



luminosity, respectively. Ilmanishi fc Wadal (120041 ) shows the positive correlation between the 



3.3 /mi luminosity and the nuclear N-band luminosity. We also obtained the same results 
for the hard X-ray luminosity. Thus, the correlation between the luminosity of the nuclear 
starbursts and central AGN is also statistically confirmed in Seyfert galaxies in hard X-ray 
luminosity, which is directly involved in AGN activity. These nuclear N-band and hard X-ray 
results definitely show that nuclear starburst luminosity is positively correlated with AGN 
power. 

For hard X-ray luminosity - 3.3 /im PAH luminosity correlation, the correlation proba- 
bility here is lower than that of the nuclear N-band luminosity. This difference may result 
from a large number of upper/lower limit data and few sample data. To analyze hard X-ray 
data in detail, we may need to eliminate these uncertainties and collect more data. We 
therefore focus on only the nuclear N-band luminosity for the following results. 

In Figure [3, we plotted the nuclear N-band luminosity versus the energy for 6.2 fim, 
7.7 /im, and 11.3 /im PAH emission. Here, we considered the 6.2 /im, 7.7 /im, and 11.3 /im 
PAH emission and plotted their energies. Since these PAH emission cover the whole host 
galaxy, we can consider them as indicators of circumnuclear starbursts. We also plotted the 
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energy for the 3.3 /im PAH emission in Fig. [7] and applied the generalized Kendall rank 
correlation statistics to these 6.2 /j,m, 7.7 /im, and 11.3 /im PAH emission simultaneously. 
The probability that a correlation is not present was found to be 0.35 % (2.9 a) for these 
PAH emission. Thus, not only nuclear starburst luminosity but also circumnuclear starburst 
luminosity are connected to AGN luminosity. 

4.4. AGN Mass Accretion and Star-Forming Activity 

In Figure El we also plotted the ratio of the nuclear N-band luminosity to the AGN 
Eddington luminosity, L N1 » 5 /L Edd versus the energy of the nuclear /circumnuclear starbursts. 
The probability that a correlation is not present was found to be 4.8 % (2.0 a) for the nuclear 
starbursts (3.3 /im PAH emission) and 16.7 % (1.4 a) for the circumnuclear starbursts (the 
other PAH emission). Thus, we statistically confirmed that the nuclear starburst is also 
positively correlated with the AGN luminosity normalized with the Eddington luminosity, 
whereas the circumnuclear starburst is only weakly correlated. 

Here PAH emission are associated with star formation, and so the horizontal axis is also 
interpreted as an indicator of star-forming activity. The vertical axis also represents the mass 
accretion rate normalized with the Eddington mass accretion rate. This ratio corresponds 
to the efficiency of gas accretion onto a given SMBH mass. Thus, these results could be 
interpreted as the relationship between the efficiency of gas accretion onto a given SMBH 
mass and the star-forming activity around the AGN. Our results show a close correlation 
between the star-forming activity that is nearer the AGN and the efficiency of gas accretion. 
Therefore, starbursts near the AGN could affect AGN mass accretion more effectively. 

5. Discussion: Link between AGN Mass Accretion and 
Nuclear/Circumnuclear Starbursts 

We found that nuclear and circumnuclear starburst luminosities are positively correlated 
with AGN luminosity. However, nuclear starburst luminosity compared to circumnuclear 
luminosity is correlated with AGN mass accretion more strongly. These results imply that the 
AGN mass accretion could be connected the star formation activity around AGN. Therefore, 
here, we discuss the interpretation of our results focusing on the effect of star formation 
activities. 



Umemura et al.l (119971 . Il998l ) showed the mass accretion onto the galactic center by the 
radiation effect from a circumnuclear starburst. They assumed a circumnuclear starburst 
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ring and a gas disk within starburst ring. 



Through radial radiation pressure from the circumnuclear starburst, the optically thin 
surface layer of the gas disk (or the whole of the optically thin disk) is made to contract. 
Then the surface layer sheds angular momentum due to radiation drag. Also, since the 
radiation drag timescale for the dusty gas has been estimated a few 10 6 yr and the radial 
radiation pressure can influence at shorter timescale further (Umemura et al. 1998), due to 
these mechanisms, luminous starbursts could carry a larger am ount of gas toward the inner 
region within the duration of the starburst phase; 10 7 yr (e.g., lEfstathiou et al.ll2000l ) and 
the typical age for AGNs; 10 8 yr (basically the Eddington timescale). 

Radial radiation pressure ma y induce star formation as well as contraction of the gas 
disk shock (lUmemura et al.lll998l ). The superbubble driven by sequential explosions of su- 
pernovae in an OB association located in the plane-stratified gas distribution may also com- 



press the surrounding g as due to shock ( Shapiro fc Field 1976 ; Tomisaka fc Ikeuchi 19861 : 



Norman fc Ikeuchilll989l ) and new stars may be born there. In fact, iMatsushita et al.l (120051 ) 
showed the clear evidence of such a self-induced starburst at the inner edge of the expanding 
molecular superbubble in M82. This effect is even more expected in the case of a luminous 
starburst. If star formation is initiated in the contracted gas, it may induce further radiation 
effects and/or a superbubble in the inner gas disk. If the starbursts occur in such a sequen- 
tial manner, the gas is carried farther into the inner region, and sequential starbursts may 
produce nuclear starbursts with some time lag (e.g., the timescale of the radiation pressure 
and radiation drag and/or the expansion time of the superbubble). 



Also, iFukuda et al.l (120001 ) showed that the gas ring in a barred galaxy could be unstable 
to gravitational instability and fragment into gas clumps. These clumps fall to the center 
by the torq ue from the m assive clumps and/or by the dynamical friction from the stellar 
component (lMakinolll9971 ). Such dynamical effects could be also related to the gas accretion 
toward the center and produce nuclear starbursts. 

Alt hough the nuclear starb urst is hidden by dusty gas, the structure is expected to be 
clumpy (IWada &: Norman! 120021 ) . In this case, not only the chaotic disturbance of the gas 
but also the radiation drag fro m the nuclear starburst may b e able to reduce the angular 
momentum of the clumpy gas (IKawakatu fc Umemura! 120021 ). Therefore, this innermost 
starburst could be directly related to the AGN mass accretion. 
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6. CONCLUSION 

To determine whether the starbursts influence AGN activities, we investigated the cor- 
relation between nuclear/circumnuclear starbursts and AGN activities. We found that cir- 
cumnuclear starburst luminosity as well as the nuclear starburst luminosity is positively 
correlated with AGN luminosity. Moreover, nuclear starburst is also positively correlated 
with the AGN luminosity normalized with the AGN Eddington luminosity, whereas the cir- 
cumnuclear starburst is only weakly correlated. This implies that close connections exist 
between starbursts nearer the AGN and AGN mass accretion. 



We thank K. Wada and T. Nagao for providing valuable discussions. We also thank the 
anonymous referee for valuable comments. YW. is supported by the Research Fellowship 
of the JSPS for Young Scientists. M.I. is supported by Grants-in-Aid for Scientific Research 
(19740109). 



Appendix 

The infrared 3- 4 /im spectra of the eight Seyfert galaxies were taken using IRTF SpeX 



(IRayner et al.l 120031 ) . The 1.9-4.2 /im cross-dispersed mode with a Of! 8 wide slit was em- 
ployed. This mode enables L- (2.8-4.1 /im) and i^-band (2-2.5 /im) spectra to be obtained 
simultaneously, with a spectral resolution of R ~ 1000. The sky conditions were photometric 
throughout the observations, and the seeing at K was measured to be in the range / . / 35-0 / . / 45 
FWHM. A standard telescope nodding technique (ABBA pattern) with a throw of 7'/5 was 
employed along the slit. The telescope tracking was monitored with the infrared slit-viewer 
of SpeX. Each exposure was 15 sec, and 2 coadds were made at each position. 

F- or G-type main sequence stars (Table [1]) were observed as standard stars, with mean 
airmass difference of <0.1 to the individual Seyfert nuclei, to correct for the transmission 
of the Earth's atmosphere and to provide flux calibration. The L-band magnitudes of the 
standard stars were estimated from their \^-band (A = 0.6 /im) magnitu des, adopting th e 



V — L colors appropriate to the stellar types of individual standard stars (ITokunagal 120001 ) . 

Standard data analysis procedures were employed, using IRAF0. Initially, frames taken 
with an A (or B) beam were subtracted from frames subsequently taken with a B (or A) 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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beam, and the resulting subtracted frames were added and divided by a spectroscopic flat 
image. Then, bad pixels and pixels hit by cosmic rays were replaced with the interpolated 
values of the surrounding pixels. Finally the spectra of Seyfert nuclei and standard stars were 
extracted, by integrating signals over 1V8-2V2, depending on actual signal profiles. Wave- 
length calibration was performed taking into account the wavelength-dependent transmission 
of the Earth's atmosphere. The spectra of Seyfert nuclei were divided by the observed spec- 
tra of standard stars, multiplied by the spectra of blackbodies with temperatures appropriate 
to individual standard stars (Tabled]). Flux calibration was done based on signals of Seyfert 
galaxies and standard stars detected inside our slit spectra. To reduce the scatter of data 
points, appropriate spectral binning is employed, depending on the continuum flux levels of 
final spectra. 
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Table 1. Observing log 



Object 


Date 


Telescope 


Integration 


P.A. 




Standard Stars 




Remark 




(UT) 


Instrument 


(Min) 


(°) 


Name 


L-mag 


Type 


Toff (K) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


NGC 931 (Mrk 1040) 


2007 August 26 


IRTF SpeX 


96 





HR 720 


4.4 


GOV 


5930 


12/im Sy 1 


F03450+0055 


2007 August 28 


IRTF SpeX 


112 





HR 962 


3.7 


F8V 


6000 


12/im Sy 1 


NGC 262 


2007 August 31 


IRTF SpeX 


80 





HR 410 


5.0 


F7V 


6240 


12/im Sy 2 


NGC 513 


2007 August 27 


IRTF SpeX 


120 





HR 410 


5.0 


F7V 


6240 


12/im Sy 2 


MCG- 2-8-39 


2007 August 30 


IRTF SpeX 


136 





HR 784 


4.5 


F6V 


6400 


12/xm Sy 2 


MCG-3-58-7 


2007 August 30 


IRTF SpeX 


80 





HR 8457 


4.8 


F6V 


6400 


12/im Sy 2 


Mrk 993 


2007 August 31 


IRTF SpeX 


112 





HR 410 


5.0 


F7V 


6240 


CfA Sy 2 


0152+06 (UGC 1395) 


2007 August 29 


IRTF SpeX 


120 





HR 508 


4.7 


G3V 


5800 


CfA Sy 2 



Note. — Column (1): Object. Col. (2): Observing date in UT. Col. (3): Telescope and instrument. Col. (4): Net on-source integration 
time in min. Col. (5): Position angle of the slit. 0° corresponds to the north-south direction. Col. (6): Standard star name. Col. (7): 
Adopted L-band magnitude. Col. (8): Stellar spectral type. Col. (9): Effective temperature. Col. (10): Seyfert 1 or 2 galaxy in the 12 (im 
or CfA sample. 
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Table 2. Properties of the Nuclear 3.3 fim PAH Emission Feature 



Object 


/nuclear-3.3PAH 


^nuclear-3.3PAH 


rest EW 3 . 3 pah 




(x 10" 14 ergs s" 1 cm- 2 ) 


(x 10 40 ergs s- 1 ) 


(nm) 


(1) 


(2) 


(3) 


(4) 


NGC 931 


<6.5 


<2.9 


<2.5 


F03450+0055 


1.8 


3.1 


1.0 


NGC 262 


<3.3 


<1.3 


<1.6 


NGC 513 


0.5 


0.3 


3.5 


MCG- 2-8-39 


<1.3 


<2.0 


<13 


MCG-3-58-7 


1.8 


3.3 


1.1 


Mrk 993 


<2.2 


<0.9 


<16 


0152+06 


<1.6 


<0.8 


<13 



Note. — Column (1): Object. Col. (2): Observed nuclear 3.3 fim PAH flux. Col. (3): 
Observed nuclear 3.3 fim PAH luminosity. Col. (4): Rest frame equivalent width of the 
3.3 (im PAH emission. 
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Table 3. Seyfert galaxies data 



Object 


z 


log ^3.3 


log £6.2 




log £7.7 


log £ll.3 


logL N1 „ 5 


log M BH 


log En 


log L x 


references 






(£©) 


(£q) 




(£©) 


(£©) 


(i Q ) 


(M ) 




(Lq) 




(i) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


Scyfcrt 1 (CfA Sample) 


Mrk 335 


0.025 


6.85 


? 




? 


? 


10.06 


6.89 


-1.34 


9.43 


1 


NGC 863 (Mrk 590) 


0.027 


<6.88 












7.58 




9.95 


2 


NGC 3786 (Mrk 744) 


0.009 


<6.36 












7.53 








NGC 4235 


0.008 


6.00 
















8.58 


3 


NGC 4253 (Mrk 766) 


0.013 


6.67 










9.71 


6.54 


-1.34 


9.09 


1 


NGC 5548 


0.017 


<6.76 










10.01 


7.92 


-2.42 


9.77 


1 


Mrk 817 


0.031 


7.36 


8.41 




8.7 


8.2 


10.43 


7.72 


-1.80 






NGC 7469 


0.016 


<6.9 


8.93 




9.27 


8.88 


9.99 


7.52 


-2.04 


9.61 


1 


Mrk 530 (NGC 7603) 


0.029 


7.24 


8.56 




8.64 


8.48 


10.01 










Seyfert 1 (12 fim Sample) 


NGC 931 (Mrk 1040) 


0.016 


<6.88 T 


<7.47 




<8.22 


7.87 


9.82 


7.13 


-1.82 


8.76 


1 


F03450+0055 


0.031 


6.91 + 


? 




? 


? 


10.07 










3C 120 


0.033 


7.4 










10.17 


7.41 


-1.76 


10.38 


2 


Mrk 618 


0.035 


7.21 










<9.95 










MCG -5-13-17 


0.013 


<6.56 










9.46 










Mrk 79 


0.022 


<7.03 


? 




? 


? 


10.15 


7.68 


-2.05 






Mrk 704 


0.030 


<6.95 










10.52 










NGC 2992 


0.008 


<6.44 










9.10 






9.30 


2 


Mrk 1239 


0.019 


<6.31 










10.11 


6.62 


-1.02 






MCG -2-33-34 


0.014 


6.53 










<8.86 










IC 4329A 


0.016 


<7.25 


? 




? 


? 


10.46 


6.64 


-0.69 


9.95 


1 


Mrk 509 


0.036 


7.36 










10.48 


7.83 


-1.86 






Scyfcrt 2 (CfA Sample) 


NGC 4388 


0.008 


<6.03 


7.68 




8.17 


7.73 


9.24 


7.22 


-2.49 


8.08 


3 


NGC 5252 


0.023 


<7.11 
















9.23 


4 


NGC 5256 (Mrk 266SW) 


0.028 


7.48 


8.96 




9.30 


8.84 


9.83 


>6.92 


< -1.60 






NGC 5347 


0.008 


<6.26 


<7.05 




<7.45 


<7.02 


9.22 


6.79 


-2.08 






NGC 5929 


0.008 


<5.59 


<6.78 




<7.32 


<6.64 


<8.37 


7.25 


< -3.39 






NGC 7674 


0.029 


7.49 


8.30 




9.03 


8.58 


10.41 


7.56 


-1.66 






Seyfert 2 (12 ^m Sample) 


Mrk 938 


0.019 


7.82 










9.85 










NGC 262 (Mrk 348) 


0.015 


<6.53 t 


<7.11 




? 


<7.29 


9.52 


7.21 


-2.20 


8.91 


5 


NGC 513 


0.020 


5.89+ 


8.45 




8.73 


8.43 


<9.69 


7.65 


< -2.47 






NGC 1125 


0.011 


6.74 










<8.95 










NGC 1241 


0.014 


<5.97 


<7.53 




<7.63 


<7.23 


<9.38 


7.46 


< -2.59 






NGC 1320 (Mrk 607) 


0.010 


6.44 










9.45 


7.18 


-2.24 






F04385-0820 


0.015 


6.61 










9.64 










NGC 1667 


0.015 


<6.21 


8.7 




9.05 


8.65 


<9.46 


7.88 


< -2.93 






NGC 3660 


0.012 


6.4 










<9.29 










NGC 4968 


0.010 


6.59 










9.49 










MCG -3-34-64 


0.017 


<6.61 










10.32 










NGC 5135 


0.014 


7.12 










9.46 


7.35 


-2.40 






MCG -2-40-4 (NGC 5995) 


0.024 


<7.26 










10.19 






9.88 


5 


F15480-0344 


0.030 


<7.13 










10.04 










NGC 7172 


0.009 


<6.28 










9.04 


7.67 


-3.14 


9.29 


2 


MCG -3-58-7 


0.032 


<7.25 










10.40 










Note. — Column (1): Object. Col. 


(2): Redshift. Col. 


(3): 


PAH luminosity at 3.3 


ftm. + Taken 


from this 


paper. Col. (7) 




7. ill. 3> and 


£jyi" 5 arc the PAH luminosity at 6.2 ^ 


m, 7.7 /^m, 


11.2 (im, 


and 


the nuclear 


N-band (10.5/xm), respectively. Col. 


(8): Black hole 


mass. Col. 


(9): Ratio of 



nuclear N-bancl luminosity to AGN Eddington luminosity. Col. (10) and (11): Absorption-corrected hard (2-10 kcV) X-ray luminosity and corresponding 
references; l. lO'Neill et al.| j2005l ); 2. iShinozaki et al.| j 2 0061 ) ; 3. |Panessa et all i2006l) ; 4. |Turner et al.| j 19971 ); 5. |Panessa fe Bass an! j2002|) . " — " shows no 
information. "?" means undetected PAH emission. 
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Table 4. Infrared luminosity 



Object 


log £lR,3.3 


l°g L ms.2 


log ilR,11.3 


£'IR,6.2/£'IR,3.3 


£lR,11.3/£ra,3.3 






(Lq) 


(Lq) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mrk 817 


10.36 


10.88 


11.05 


3.30 


4.94 


NGC 7469 


< 9.90 


11.40 


11.73 


> 31.52 


> 68.21 


Mrk 530 (NGC 7603) 


10.24 


11.03 


11.33 


6.14 


12.41 


NGC 931 (Mrk 1040) 


< 9.88 


< 9.94 


10.72 




> 6.98 


NGC 4388 


< 9.03 


10.15 


10.58 


> 13.14 


> 35.80 


NGC 5256 (Mrk 266SW) 


10.48 


11.43 


11.69 


8.88 


16.36 


NGC 5347 


< 9.26 


< 9.52 


< 9.87 






NGC 5929 


< 8.59 


< 9.25 


< 9.49 






NGC 7674 


10.49 


10.77 


11.43 


1.90 


8.79 


NGC 262 (Mrk 348) 


< 9.53 


< 9.58 


< 10.14 






NGC 513 


8.89 


10.92 


11.28 


106.79 


247.67 


NGC 1241 


< 8.97 


< 10.00 


< 10.08 






NGC 1667 


< 9.21 


11.17 


11.50 


> 90.89 


> 196.73 



Note. — Column (1): Object. C olumn (2): Infrared luminosity of the nuclear starb urst. LrR.3.3 = L3.3 X 10 3 
iMouri et al"lll99d ; Ihnanishill2002h. Col. (31 and (4): Ltr, 6 .2 = £-6.2 x 3 X 10 2 ) llPeeters et al.l l2004h and 
^IR 11.3 = ^11. 3 X 7 X 10 2 ISoifer*^t al.l20o3) are the infrared luminosity of total starbursts, which were estimated 
from 6.2 fim and 11.3 fim PAH luminosity, respectivery. Col. (5) and (6): Infrared luminosity Ratios of total 
starburst to nuclear starburst. " — " shows that both the infrared luminosity of the nuclear starburst and the 
total starburst are upper limit data. 
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Fig. 1. — The schematic view of the nuclear /circumnuclear starbursts. The nuclear starburst 
(within a few 100 pc) and the circumnuclear starburst (several 100 pc - a few kpc) are 
distributed around AGN. 
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Fig. 2.— Zoom-in spectra around the redshifted 3.3 /zm PAH emission feature (lower 
arrows) of the observed eight Seyfert galaxies. The horizontal and vertical axes are the 
observed wavelength in /iin and F A in 10~ 15 W m -2 /im -1 , respectively. 
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Fig. 3. — 5.2 — 14.5/iin spectra of Seyfert 1 galaxies. The horizontal and vertical axes are 
the observed wavelength in /zm and F\ in Jy, respectively. The solid straight lines are the 
adopted continuum to estimate each PAH emission fluxes. The solid curves are Gaussian 
profile fittings of PAH emission. The lower arrows and the question marks indicate the weak 
(lower than 3a) and undetected PAH emission, respectively. 
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Fie;. 4. — Same as Fie;. [3l but these spectra are Seyfert 2 galaxies 
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Fig. 5. — The horizontal and the vertical axes are the nuclear N-band luminosity and the 
absorption-corrected hard (2-10 keV) X-ray luminosity, respectively. The percentage repre- 
sents the probability that a correlation is not present. 



-26- 



11 

10.5 
10 
■ 9.5 



J 9 

60 



7.5 < 

11 
10.5 
1 10 

j 9 

CO 

■S 8.5 
8 
7.5 




5.5 



5.5 



6.5 7 

l°gL 3 .3p AH [L s „ 



7.5 



•A 
A 



5.3% Se y fert 1 * 

Seyfert 2 * 



6 6.5 7 



7.5 



8.5 



Fig. 6. — upper. The horizontal axi s is nu clear 3.3 /im PAH emission luminosity detected 
inside slit spectra (jlmanishi fc Wadal (120041 ) and this paper). The vertical axis represents the 
nuclear N-band (10.5 um) luminosity mea sured with ground-based two-dimensional camera 
with a 1".5 aperture (IGorjian et al.ll2004l ). bottom: Same as upper figure, but the vertical 
axis represents the absorption-corrected hard (2-10 keV) X-ray luminosity. The percentage 
represents the probability that a correlation is not present. 
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Fig. 7. — The horizontal axis is the energy of the 3.3 /zm, 6.2 /im, 7.7 /im, and 11.3 /im PAH 
emission and the vertical axis is the nuclear N-band luminosity. The percentages represent 
the probability that a correlation is not present. 
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Fig. 8. — Same as Fig. but the vertical axis represents the ratio of the nuclear N-band 
luminosity to the AGN Eddington luminosity. 



